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Abstract 

In the present article we perform analytical and numerical calculations related 
to persistent currents in 2D isolated mesoscopic annular cavities threaded by 
a magnetic flux. The system considered has a high number of open channels 
and therefore the single particle spectrum exhibits many level crossings as the 
flux varies. We determine the effect of the distribution of level crossings in 
the typical persistent current. 

INTRODUCTION 

Recent developments in micrometer-scale technology have made possible the fabrication 
of devices small enough to make electrons phase-coherent inside. Therefore the transport 
properties of those samples at very low temperatures exhibit features characteristic of quan- 
tum coherence of the electronic wave function along the sample. The observed Aharonov- 
Bohm oscillations in the resistance of a loop pierced by a magnetic flux is one of the most 
relevant manifestations of the interference phenomena typical of phase coherence [1].[2|] . 
Although the transport properties have been intensively investigated both theoretical and 
experimentally, the persistent current problem remains less understood. In 1983 Biittiker, 
Imry and Landauer, following earlier work by Byers and Yang on superconducting rings || , 



proposed the existence of such currents in one dimensional mesoscopic loops in the presence 
of weak magnetic fields 0] . These equilibrium currents are a consequence of the nature of 
the eigenfunctions and their flux sensitivity, which is strictly of the Aharonov-Bohm type. 
The current is a periodic function of the magnetic flux with fundamental period <po = hc/e. 
In ideal clean rings at K (ballistic regime) the current is proportional to Iq = e vpjL 
where vf is the Fermi velocity and L is the length of the loop. On the other hand, in the 
diffusive regime (when impurities are present and the elastic mean free path / is shorter than 
the typical sample size) theoretical studies predict that the current decreases by a factor l/L. 
The problem of persistent currents did not deserve much attention until recent years when 
experimental measurements have been performed on metallic rings in the diffusive regime 
|5||| and in a few-channel mesoscopic semiconductor ring 0. For the metallic systems, the 
experimental value of the current is ~ (0.3 — 2) Jo, 1 or 2 orders of magnitude greater than 
the theoretical predictions consistent with diffusive motion, that is with no apparent correc- 
tions due to disorder. The value of the current reported by Mailly and collaborators on the 
mesoscopic ring, is also of the order of Iq, but in their experiment l/L ~ 1. 
Owing to the discrepancies between the observed and predicted values of the current, much 
of the theoretical efforts were put to find a mechanism which could compensate the effect of 
impurities. We are not going to describe those approachs here. For a recent review of the 
problem, see |§. 

Most of the theoretical studies mentioned above where restricted to one dimensional geome- 
tries ||. The effect of the number of channels in the persistent current, although investigated 
by many authors, is not yet completely understood. The computations of these currents in 
2D (multichannel) geometries is much more difficult, and all the studies for multichannel ge- 
ometries have been performed employing discrete models or cylindrical geometries in which 
the transverse channels are decoupled for the conduction mode [|T0|-|T2"fl . 

In this paper we deal with the persistent current problem for a system of N non- 
interacting electrons confined in a mesoscopic clean annular 2D geometry at K. This 
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geometry is useful to describe the real metallic "rings" employed in most of the experiments 
mentioned above. As an example in Ref. J7[ the ring has internal diameter 2/im and external 
diameter 3.4/im, conforming a 2D cavity. 

Our system is formally a quantum billiard whose hamiltonian H depends on a parameter 
(the magnetic flux). The variation of the parameter preserves the commutator [H, L] = 
(L being the angular momentum), and therefore the single particle spectrum displays large 



amount of crossings as the parameter is varied |13[ . Such degeneracies and their features are 
fundamental ingredients that affect the persistent currents when many channels are open. 
In fact, they are not only responsible for the multiple jumps that the persistent currents 
show as a function of flux (for fixed number of particles N) but also for the large oscillations 
in the typical current as a function of N after averaging on the flux. These oscillations 
(which are inherent to the system because they depend on the features of the crossings of 
the single particle energy levels) make difficult to establish an average behavior with N. 
In the present work we show that not only the actual average behavior but also the large 
oscillations emerge when the statistical properties of the degeneracies are considered. 
The paper is organized a follows. In Sec. |, we present our system (we call it the Aharonov- 
Bohm annular billiard). Sec. [TT] is devoted to summarize some results and properties of 
persistent currents in ID and 2D. In Sec. |TTT| , we introduce the properties of level crossings 
that allow us to perform a statistical approach to the problem. Section |V| is addressed 
to exploit the properties of the level crossings in order to determine the dependence of the 
average behavior and the fluctuation of the persistent current on the number of particles N. 
Finally, in Sec. [V] we present the concluding remarks. 

I. THE AHARONOV-BOHM ANNULAR BILLIARD 

We consider a 2D structure with the geometry of an annulus. See Fig. [I]. The cavity 
consists on the planar region limited by two concentric circles of outer and inner radii R and 
r respectively. In the following we will take the area equal to 7r and we define the parameter 
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A = R/r. The spatial degrees of freedom are the azimutal angle 9 and the radial coordinate 
p which varies between r and R. 

Let N be the number of non-interacting electrons on the cavity. We are interested in the 
persistent current caused by the application of an homogeneous time independent magnetic 
flux threading the cavity axially. We disregard the effects of the magnetic field piercing 
the body of the annulus, so the electrons move with uniform rectilinear motion inside the 
cavity. We choose the gauge in which A = (j)/{2irp) 9 where 9 is the azimutal unit vector. 
The single particle spectrum results from the eigenvalue equation 

A* + — -—-—* + k 2 V = 1.1 

p z o9 p z 

where A is the Laplacian in polar coordinates. We define the scaled flux a = 4>/(po and we 
use units such h 2 /2m = 1, so the energy is E = k 2 . 

We apply Dirichlet boundary conditions at p = r and p = R and periodic boundary 
condition in the azimutal direction. 

The Eq. ( |1 . 1| ) is separable in polar coordinates and we factorize ^(p, 9) = T{p) expimO 
with m = 0, ±1,±2... the orbital quantum number. The eigenfrequencies k Ujn results from 
the solution of the equation 

J v {z\)N v (z) - J u (z)N u (z\) = , 

v = m + a , (1.2) 

where we have defined z = kr and n = 1,2..., is the radial quantum number. J u and 
N v are the Bessel functions of the first and second kind, respectively. The corresponding 
eigenfunctions *$f(p,9) are : 

V Vjn (p,6) = A Vtn [J v (k Vtn p)N v (k v>n r) - J v {K,nX)N v {k^ n p)} exp(im9) (1.3) 

where A UjU is the normalization constant. 

All the eigenstates and all the equilibrium physical properties of the system are periodic 
functions of the flux with period 0o @- Moreover as the energy spectrum is symmetric with 



respect to = 4>o/2, in the following the parameter a will take values between and 1/2. 
For the states with m and — m are, in general, not degenerate. This is the reason for 

the existence of an equilibrium persistent current. 

We stress that for the present system E Ujn can not be written down as a simple function 
of the numbers v and n as it happens, for example, in the cylindrical geometries in which 
the eigenenergies are cuadratic functions of the quantum numbers. This constitutes the 
fundamental difference between the annular cavity and other 2D geometries studied so far, 

H. 

To obtain the eigenenergies E v>n we have numerically solved Eq. ( |1.2|) . For each fixed 
value of A, m and n, we have taken six equally spaced values of a between and 1/2. Then 
we have verified that the best quadratic fit was quite satisfactory for the considered values 
of A. Fig. |2| shows a region of the energy spectrum as a function of a for A = 10 obtained 
by the described procedure. In the following we will consider 

E m + a ,n = A min (\)a 2 + B m>n (X)a + C m>n (X) . (1.4) 

Let us remark that for the cylindrical geometries (ie. a cylindrical surface of area L x 
L y ) the coefficients can be determined exactly and they are, A = (2ix / L) 2 (without any 
dependence on the quantum numbers) and B = 2 m (2tt/L) 2 which depends only on the 
orbital quantum number m. The other quantum number appears only in the constant term 
C. 

II. PERSISTENT CURRENTS IN N ELECTRON SYSTEMS 

In the present section we summarize the main results concerning persistent currents in 
ID (rings) and 2D (cylinders and annula) systems. 

For a system with iV non-interacting electrons, the total current is It = J2n=i In @- 

In order to characterize the typical current we define I typ = J (I T ) a , where the symbol (..) a 

means flux average. This is 
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I% v = 2 fj 2 4 da ; (2.1) 
where we have profited the symmetry of the spectrum with respect to a = |. 

A. ID Systems 

For a ID ring of circumference L threaded by a magnetic flux 0, the velocity of an 
electron in the state with orbital quantum number m is v m = jdE m jdk m . Therefore as 
k m = 2tt (m + a)/L, the current carried by the level E m results ||] 

uEm 

/„ = ~ ■ (2.2) 

It is well known that for the ID ring geometries, the crossings between levels in the 
range < a < 0.5 occur only for a = or a = 0.5. This fact is a direct consequence of the 
existence of only one channel in the system. The total current for the N-electron system is, 

N a N odd 

(2.3) 

j§ti N (2 a - 1) N even 

where we have explicitely used kp = tcN/L. Therefore the typical current is I typ ~ N. This 
is consistent with the well known result for ID geometries, It yp = Io with Iq = e Vp/L 0. 



B. 2D Systems 

In the case of 2D cylindric geometries, as the energy of an individual state is separable 
in two terms, one associated with the conduction mode and the other with the transversal 
one, it is straigthforward to verify Eq. ( |2.2j ). For the annular geometry, that equation is not 
obvious. In order to prove Eq.( p.2|) we begin computing the current density j : 

j = -— [** (-i h V + -A) * + * (i W + -A) **] . (2.4) 

From the Eq. (|1 . 1|) it is straightforward to show that 
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dH e 



dip m c 



[-1 hV + -A\. 



dl 

~d0 



(2.5) 



Comparing Eq. (|2.5| ) and Eq. (|2.4| ), we obtain 

2e -t dl T ^dH 

— 3 ■ = * 

c d(j) d(f) 




Taking into account the functional form of the vector potential A and the Hellmann-Feynman 



theorem |14], the current carried by the state i results: 



h = - e j e dp= -c 



(2.6) 



d<f) 



where Ei = k/ 2 and the limits of integration are r and R respectively. 

Therefore, following Eq.( |2.6| ), the total persistent current It, is the current through a tran- 
verse section of the annulus. 

For 2D geometries, as we have already mentioned, the effect of the number of channels 
on the typical currents is not yet fully understood. The number of channels M is defined as 
the maximun value of the transverse quantum number inside the Fermi surface at zero flux. 
We define the value of Fermi energy at zero flux as Ep = En, En being the energy of the 
highest ocuppied level. We emphasize that in the considered isolated cavities the number 
of electrons iV is constant for all values of the flux. For values of a ^ 0, the average Fermi 
energy is defined by Weyl's formula (E F ) = N/4 (in the stipulated units) 
One important point is to know how M depends on the number of particles iV and on the 
parameter A. Employing the expansion of the zeros of the radial functions (solutions of the 



where the symbol [..] means integer part. 

For A — > 1, the second and third terms in the last equation go to and we obtain the trivial 
result for ID geometries, M = 1 independently of the number of particles N. 

Most of the experimental relevant situations involve many particles. For values of A > 1 
and for N » 1 the second term in the last expression dominates and M reads: 



Eq.flig)) we obtain [|TJJ: 




7 



M 



2 V^V (A-l) 



(2.7) 



Fig. |3] shows the total current vs flux for N = 2747 and A = 10. This corresponds to 
M = 60. The existence of many channels implies that, opposite to what happens for ID 
geometries, crossings between levels of different channels occur for < a < 0.5. 



where the coefficients A's > 0, and the coefficients B's have alternated signs. The sum 
extends to all the filled levels N. The total current has two contributions, the first term 
which has an explicit linear dependence on a and the second one without explicit dependence 
on the flux. Nevertheless, Fig. |^ shows the total current which looks like a sawtooth. Between 
two succesive crossings the curve is a linear function with the same positive slope in each 
interval. The negative jumps in It are for values of the flux a[~^ where crossings between 
states occurs. This qualitative behavior can be understood as follows. Let us suppose that 
a given state is occupied before a crossing and empty after it. This corresponds to replace 
in the Eq.fl2.8p one value of A m>n by another one. As iV >> 1, the last replacement gives 
essentially the same slope. On the other hand, the second term is strongly affected because 
a positive value of B nm is replaced by a negative one, and this produces a significative 
negative change in the value of the non homogeneous term of It- Let us remark that 
between two adjacent a\ ^ 's there is another crossing at a[ + ^ that is irrelevant for It- We 
denote n^\n^) the total number of crossings at the positions As it will be 

clarified in the following sections, the features of the crossings between levels are strongly 
related to the shape and the magnitude of the total and typical currents. 



Employing the expansion given by Eq. ( |1.4p , the total current for the N electrons system 



can be written as : 




(2.8) 
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III. THE DISTRIBUTION OF CROSSINGS 



We devote this section to establish some relevant properties of the crossings between 
states for integrable hamiltonian systems that depend on a single parameter. Following the 
semiclassical expansion for the density of states developed by Berry and Tabor for integrable 



systems ||17|| , we have started writting the density of crossings per unit of energy and flux 



(for a given value of the parameter A) as 

p c (E, a, A) = p c + p, 



osc . 

C J 



p c is the average density of crossings (where the average is taken over energy and flux) and 
p c osc is the fluctuating part. We have demonstrated that 

p c = C(A)v^, (3.1) 



and is independent of a [18|]. Therefore the average number of crossings that a given level 



n has between < a < 0.5 is h c = 0.5 p c . Moreover, 

n (+) = n(") = n c /2 . (3.2) 

Fig. ^ shows n c vs. n, n c being the number of crossings of the level n, for A = 10 and 
N = 3650, together with the ensemble average n c taken on a range 1 << An < N. The best 
power fit to n c is n c = 0.6182 n 13 with j3 = 0.50441, in strong agreement with our theoretical 
prediction. It is interesting to note that while the number of channels scales as M ~ \fN 
(see Eq.( [2.7|) ), the mean number of crossings for the highest occupied level is proportional 
to M . That means that on average, the last occupied level experiments as many crossings 
as there are channels present in the system. The higher fluctuations observed in the Fig. |] 
correspond to the apparition of a state with m = (at a = 0) in the energy region consid- 
ered. 

Another important result concerns the way in which the crossings for a given level n are 
distributed in the interval < a < 0.5. We define the normalized nearest neighbour spacing 
as Si = fi c {a.i — ai-i), where «j denotes the position of the crossing i. Fig. || shows the 

9 



nearest neighbour spacing distribution P(s), obtained after an exhaustive numerical com- 
putation. This results in a Poisson distribution as can be checked in the mentioned figure, 
where the numerical result is displayed along with the exact Poisson distribution. In other 
words, for a given level n with n c crossings, the corresponding a' { s can be seen as a sample 
of n c uncorrelated random variables in the interval < a < 0.5. Therefore, if we generalize 
the above definitions to n-nearest neighbour spacings as Si = h c (ai+ n — the associated 



spacing distributions are |19 



P(n, s) = — j exp (— s) . (3.3) 

In the notation of the last equation P(s) = P(0, s). Such distributions of the spacings de- 
termine the fluctuations of the typical current, as we will show in the forthcoming section. 



IV. TYPICAL CURRENTS OF 2D SYSTEMS AND THE CROSSINGS AS 

RANDOM VARIABLES 

At this point we have all the ingredients to quantify the way in which the distribution of 
crossings determines the scaling of the typical current with the number of particles present 
in the system N. 

We want to compute if (see Eq. (|2.1|) ). Taking into account the definition of the total 



current Eq.( [2.8|) , and the comments stated at the end of Sec. [IIB| , we can write the square 
of the total current for < a < 0.5 as, 

i| = (A T a-F(a,ai- y )) 2 , (4.1) 
Ffoof-^EAieCa -a, H ), (4.2) 

i=l 



where we have defined At = 2 X^m,n (see Eq. ( |2.8|) ) and 6(x) is the step function. The 
positive definite quantity A; is the absolute value of the jump that the total current has 
at the position af' . As we mentioned before, all the jumps have the same negative sign, 
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therefore we have written down a minus sign preceding the second term of Eq. (|4.1|) . We 
stress that in the last equation denotes the total number of crossings (for < a < 0.5) 
for which It shows jumps. For the sake of symplicity in the notation, in the following we 
drop the supraindex in the a\ ^ and in n^~\ It yp , and obviously its squared value, can be 
written function of the random variables a,-, 

(A 2 nc /l \ n ° /l \\ 

+ E A * ^ (2 - Max ^ a s)) - A ? £ A < (4 -«m » ( 4 - 3 ) 

where Max (a, 6) means the maximum among a and 6. Figure ^| shows the exact numerically 
obtained I 2 as a function of the number of particles in the system N for A = 10. The 
extremely fluctuating behavior with N is displayed. 

Our aim is to determine the average dependence of It yp 2 with the total number of particles 
in the system N. This implies performing an average over N when the number of particles 
varies in a range 1 << 5N < N (N-average). In the same figure a solid line shows this 
average for w 2000. In the present case, such value of 5N is the required to obtain a 
smooth curve. 

The fluctuations in the fig. |] are due to the different patterns of the sequences of a[s that 
appear when N is varied. Therefore, we assume that the N-average is equivalent to statistical 
averages over different realizations of the random sequence of a[s. We define the average 
typical current as 

i 2 = <4p> ( 4 - 4 ) 

where the symbol (..) means average over different realizations of a^s. 
I characterizes the magnitude of the current that is relevant for the experimental investiga- 
tions 0. 

Before computing J 2 , we will give a qualitative explanation of the origin of the fluctuations 
in lf yp . The total current grows linearly between crossings and therefore its growth is in 
direct proportion to the distance between succesive crossings (see Fig. |3|). There are two 
extreme realizations of the sequence of a[s that limit the expected behavior of I 2 . 
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Let us first suppose that the crossings conform a regular sequence of equally spaced a[s 
(what is known as a picket-fence sequence) and Aj = A VI In this case it is easy to verify 
that I? ~ N/ric- As we have already stated in Sec. [II], n c ~ ^/N (in this case n c = n c ) 
and therefore Ifyp ~ 

N. 

The second realization corresponds to a sequence of a[s in which one distance between cross- 
ings is much larger than the other ones. This resembles the features of the ID systems for 
which in fact, there is only one crossing in the interval < a < 0.5. It is straigthforward to 
show that for the realization considered Ityp ^ 

N 2 . 

For the annular billiard, as we have shown in the previous section, the nearest neighbour 
spacing distribution between crossings is Poisson. Therefore one expects the actual behavior 
of I^yp to be highly fluctuating with N around a mean value ~ iV 7 with 1 < 7 < 2. At the 



end of the present section we will see that the Eq. ( f4.3|) , that depends on the particular 



realization of the sequence a[s, fullfills the above expectation. 

In order to proceed with the computation of J 2 , it is necessary to assume some behavior 
for the amplitudes Aj. In the Fig. |3|, it can be seen that the amplitudes of the jumps are 
not essentially very different on average. This fact suggests us to assume 

A, = A = ^ n c , (4.5) 

Vi. This ansatz will be strongly justified by our results. 

From Eq. (|4.3|) and performing the ensemble average over a^s we obtain: 



I 2 = 2 ( § - ^ A n c + A T A I a%_ x ) j + ^ n\ - A 2 f £(2 j - l)(a 2 ,_ 1 ) 



(4.6) 



where the sums run over the odd values because we have assumed that the odd crossings 
contribute to the jumps in the total current (it is worth to note that the election is arbitrary 
and an equivalent treatment assuming jumps in the even crossings gives the same results). 
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To evaluate (o^-i) and (a^-i) we have considered a%j-\ = s^j-i, where S2j-i is the 
spacing between the crossing in a = and the crossing in a = OLij-\- Therefore, employing 
the distributions Eq. ( |3.3j ) (adequately normalized to nearest neighbour mean spacing d = 
1/(4 n c )) we have obtained: 



J2(2j-l)(a 2j -i) = d 

3=1 



4-3 1 ~ 

- n, — - n, 
3 c 3 



here we have employed 



(a 2j -i) = d (2j - 1) , 

(a 2 2j _ 1 ) = d 2 2j(2j~l). (4.8) 



After replacing Eq. (|4.7Q in Eq. ( |4.6| ) we finally obtain : 

I = 2 \ — + ^ + — — + T — . 4.9 

\ 24 6 6 16 12 48 n c / K J 

The above expression allows us to classify each term according to the dependence on N. 
Taking into account the Eqs.( [3.1D , ( p.2|) and (|4.5|) and the fact that At oc N, we can verify 
that each of the first three terms in the Eq. ( |4.9|) are of order 0(N 2 ) and that their total 
contribution vanishes. The fourth term is of order 0(N 3 ^ 2 ) and the two last ones are O(N). 
We have obtained for the mean typical current, 

A 2 A 2 

p = Ar_ + ^L_ (4.10) 
16 n c 48 n 2 y ' 



and therefore for large N, 



0(N Z ' 2 ) . 



Eq. ( |4.10| ) constitutes one of the main results of the present article. We have plotted in 
the Fig. |6|(a), besides the exact numerical results (dots), the N-average with 5N = 2000 
(solid line) togheter with its best power fit ~ iV 7 (dashed line). The resulting value for the 
exponent is 7 = 1.498 in strong agreement with our theoretical result. 
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At this point, let us mention that the choice of SN = 2000 to perform the N-average could 
resemble somehow ambiguous. However, if we perform an N-average with 5N < 2000, 
although the resulting curve displays oscillations, the best fit leads to the same scaling 
with N (see Fig. ©(b)). It is important to remark that the Eq. (fO|) can be employed to 
evaluate the typical persistent current for a given realization of the a[s. In particular for 
the picket-fence sequence (note the absence of the ensemble average in the following), 



4 1 
Z!(2 j ~ 1) «2j-i = d y - hi - - n, 



YAj-x = # { 3^e - 3 n c ) . (4.11) 

The last equations lead to a cancellation of all the orders greater than 0(J\f) in the Eq. ( [4 . 3| ) 
for If yp . That is 

A 2 

*w 48 n\' [ ] 

Equation ( |4.12j ) satisfies I^ yp ~ in accordance with our previous qualitative discussion 
(see the paragraph below Eq. Q4.4D). 

On the other hand, we expect some other particular realizations of the sequence of a[s 
for which there is no cancellation of the terms of order 0(N 2 ) in Eq. ( ^4.3| ). In this case 

n w - n 2 . 

We stress that the considered particular realizations manifest themselves when I 2 yp takes 
the lowest and highest values respectively (see Fig. H(b)). 



V. CONCLUDING REMARKS 

In the present work we have characterized the behavior of the persistent current for a 
system of non-interacting electrons confined in a mesoscopic bidimensional clean cavity as 
a function of the number of particles in the system N. We have deduced relevant properties 
of the single particle spectrum such as the scaling of the number of open channels M and 
average number of level crossings h c with N. 
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The main goal of this article has been to determine the iV-average behavior of the typical 
current through the statistical properties of the degeneracies present in the single particle 
spectrum as a function of the flux, when many channels are relevant. We can separate the 
dependence on N of the typical current in two contributions. The first one is the average 
value taking on a range 1 << 5N < N and constitutes the smooth contribution. The second 
one, depends on the particular realization of the sequence of crossings a[s (the values of flux 
where degeneracies in the spectrum occur) and leads to an extremely fluctuating behavior 
with N. We have related this fact to the statistics of the distribution of crossings which, as 
we have determined is Poissonian. The smooth part of the typical current has been obtained 
under the assumption that the ensemble average over different realizations of the random 
sequence of crossings a[s is equivalent to the N-average. The obtained dependence on N 
of this smooth contribution is I 2 ~ N 3 ^ 2 which is in strong agreement with the numerical 
results (see Fig. |^). On the other hand we have evaluated the typical current for particular 
realizations of sequences of a[s that lead to completely different dependence on N, namely 
tfyp ~ and I 2 yp ~ N. These realizations are responsible for the large oscillations of I 2 . 
Last but not least, some words of caution are in order in connection with experimental 
observations of persistent currents. It is necessary to emphasize that, owing to the extreme 
sensitivity of the typical current on the particular realization of a-s, any experiment devoted 
to establish the scaling of the average typical current on N (or on the number of open channels 
M) in mesoscopic 2D cavities would involve measurements on the order of a thousand of 
samples (remember that SN « 2000, for the annular geometry). 
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FIGURES 




FIG. 1. Annular billiard threaded by a magnetic flux <p. The internal and external radii 
denoted r and R respectively. 
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FIG. 2. Single particle energy levels as a function of the adimensional flux a, corresponding to 
the lowest region of the spectrum for the annular billiard with A = 10. The circles correspond to 
exact values obtained from the zeros of the cross products of Bessel functions. The small dots lines 
are the quadratic fits. See the text for details. 
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FIG. 4. Number of crossings n c for < a < 1/2 as a function of the number of levels N. The 
points are the exact values numerically calculated. The thin line corresponds to an average over 
AiV = 100 while the dashed line corresponds to the best power fit (n c = 0.61822 x jV - 50441 ) to 
this average. 
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FIG. 5. (a) Nearest neighbour spacing distribution P(s) numerically calculated. We also show 
the exact Poisson distribution P(s) = exp— s. (b) Log-lin plot of P(s). 
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FIG. 6. (a) The square of the typical current lf yp as a function of N. The points correspond to 
exact values numerically calculated. The solid line corresponds to an N - average with AN = 2000, 
while the dashed line is the best power fit lf yp = 2.009 x jV L4966 . (b) The same exact results 
of a) but the iV-average was performed with AiV = 1000. In this case the best power fit is 
lf yp = 1.7405 x iV 1 - 5558 . Here, we have drawn the limit curves (oc N 2 and oc N) for particular 
realizations of sequences a[s. 
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